Solid-liquid interfacial premelting 
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We report the observation of a premelting transition at chemically sharp solid-liquid interfaces 
using molecular-dynamics simulations. The transition is observed in the solid-Al/liquid-Pb system 
and involves the formation of a liquid interfacial film of Al with a width that grows logarithmically 
as the bulk melting temperature is approached from below, consistent with current theories of pre- 
melting. The premelting behavior leads to a sharp change in the temperature dependence of the 
diffusion coefficient in the interfacial region, and could have important consequences for phenom- 
ena such as particle coalescence and shape equilibration, which are governed by interfacial kinetic 
processes. 



The term premelting refers to the formation of ther- 
modynamically stable liquid films at solid interfaces sub- 
jected to temperatures below but near the bulk melting 
temperature (T m ) 0, Q. Surface premelting (SP), the 
formation of a premelting layer at a solid-vapor inter- 
face, was first observed experimentally by Frenken and 
van der Veen Q using proton scattering. Since then, 
advances in experimental techniques have provided pow- 
erful tools for direct atomic-resolution observations of 
this surface phase transition These experimental 

studies have been complemented by a number of de- 
tailed atomistic simulations probing the existence and 
atomic-level mechanisms of SP (Refs. Illl - fl8l and refer- 
ences therein). Premelting at solid-solid interfaces has 
also been reported in the literature and can take two basic 
forms: premelting at solid-solid hctcrophasc boundaries 
and grain-boundary premelting (GBP). Examples of the 
former are found at interfaces between Pb and Alflij]. 
SiO and A1 2 03[2(J, as well as in ice at solid substrates, 
which plays a role in frost heave [III (UHl^]. Premelt- 
ing at grain boundaries has been the subject of numer- 
ous continuum [24l - [3(ij , atomistic (Ref. [H and references 
therein) and experimental 3l| studies. 



In this work, we report molecular-dynamics (MD) sim- 
ulation results on Al/Pb solid- liquid interfaces that pre- 
dict the existence of a third class of premelting transition, 
namely, solid-liquid premelting (SLP). In this process a 
premelting liquid layer forms at a solid-liquid interface, 
below T rn of the solid. Such a transition should, in princi- 
ple, be possible at chemically heterogeneous solid-liquid 
interfaces in which the melt phase of the solid and the liq- 
uid phase are mutually immiscible, as is the case with Al 
and Pb near T m for Al. Solid-liquid premelting has not 
been previously reported under ambient pressures, cither 
experimentally or by simulation - although some evidence 



of SLP have been reported in simulations of solid-liquid 
interfaces under extremely high pressure (diamond anvil) 
32 . [33J]. As we will demonstrate, the process of SLP 
leads to a sharp change in the temperature dependence 
of the diffusion coefficient in the interfacial region, which 
is expected to have important consequences for the ki- 
netics of interface-controlled processes such as particle 
shape equilibration or coalescence of liquid nanoparticlcs 
in solids governed by interface-mediated Brownian mo- 
tion jalEl- 

Thermodynamically, premelting occurs near T m when 
the interfacial free energy, 7 SQ , between the solid and 
another phase a (a — solid, liquid or gas) is larger than 
the sum of that for the solid- melt (7 s i) and melt-ev (71a) 
interfaces: 



A7 = 7 S q - [7sl + 7la] > 



(1) 



Thus, if the undercooling (AT = T m —T) is not too great, 
it is thermodynamically favorable to form a thin film of 
metastable liquid because the increase in bulk free en- 
ergy is more than compensated by a lowering of the total 
interfacial free energy. The width of the interface as a 
function of undercooling depends on A7 and AT, as well 
as the nature of the potential of interaction between the 
two interfaces - the so-called " disjoining potential" . The 
pioneering theoretical studies by Kikuchi and Cahn [36[ 
on grain-boundary premelting and Lipowsky 37( on SP 
both predict a logarithmic dependence HU of the width 
of the premelting layer with respect to the undercool- 
ing IH: 



w(AT) = -w In [AT /T ) 



(2) 



where u>o and T) are constants specific to the given in- 
terface. 
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The System: The Al-Pb system is an ideal model alloy 
for the study of chemically heterogenous solid-liquid in- 
terfaces. The phase diagram is a simple monotectic that 
has a broad liquid-liquid miscibility gap, negligible sol- 
ubility of Pb in the Al solid phase and a large melting 
point separation (600 K for Pb and 933 K for Al). We 
have previously reported results from MD simulations on 
this system at 625 K, a temperature just above the melt- 
ing point of Pb. [!(| The simulation results are consistent 
with an in situ transmission-electron-microscopy (TEM) 
study of liquid Pb inclusions embedded in a crystalline 
Al matrix |41( in that the (111) interface is shown to be 
faceted, while (110) and (100) are rough at this tempera- 
ture. The experiments also show that the (111) interface 
undergoes a roughening transition about 100 K below the 
melting point of Al. 

Simulation Details: In our simulations of the Al/Pb 
solid-liquid interface, we employ a classical many-body 
potential developed by Landa, et al. 42j to model the 



interatomic interactions. This potential predicts mutual 
immiscibility of Pb and Al in the solid state and a large- 
liquid state miscibility gap consistent with the experi- 
mental phase diagram up to 1200K. The melting points 
of this potential are 615.2(2)K and 922.4(2)K for Pb 
and Al, respectively. The MD simulations arc performed 
using LAMMPS [43[. Equilibrated SL interfaces are set 
up at varying temperatures, T, ranging from 625 to 900K, 
separated by 25K up to 90K with additional simulations 
at 912, 920, 921 and 922K. Three crystallographic ori- 
entations for the Al-Pb interfaces are examined: (100), 
(110) and (111). 

To produce the equilibrated interfaces, constant-area, 
constant- normal-pressure MD (NP Z AT) simulations up 
to 50 ns in length are used to yield the appropriate equi- 
librium number density, p, pressure (1 bar) and composi- 
tion. These are followed by constant NVT simulations to 
collect production data. Five replica systems (each con- 
taining two independent interfaces) are used at each tem- 
perature and orientation to improve statistics. For addi- 
tional details as to the methods of interface set-up, equi- 
libration and analysis, see the supplemental materials [44j 
and Ref. 0. 

Characterization: The solid-liquid interface is char- 
acterized through the determination of interfacial pro- 
files, which show the change in specific properties (e.g., 
density, local structural order, composition and diffusion 
constant) as functions of the distance normal to the inter- 
facial plane, defined as the z direction. The z-coordinatc 
is measured relative to a Gibbs Dividing Surface (GDS), 
defined here such that the excess number of Al atoms. 



ing layer, we utilize two different order-parameter (OP) 
profiles. The first profile uses a local structure OP that 
distinguishes solid from liquid phases fjij and is normal- 
ized to be 1 in the solid phase and in the liquid phase. 
The second is a compositional order parameter equal to 
1 in a pure Al system and for pure Pb. For a Al/Pb 
solid-liquid interface without a prcmclting layer, these 
two order-parameter profiles will be approximately coin- 
cident; however, in the presence of a prcmclting layer (a 
liquid Al layer separating solid Al from liquid Pb), the in- 
terfacial position indicated by these two OP profiles will 
be separated by the width of the premelting layer, w. 




Tai is zero. 44[ To determine the extent of the premelt- 



FIG. 1. Snapshots of Al-Pb (111) solid-lilquid interfaces at 
three different temperatures. Top panel: 625K, interface is 
faceted [4(J; Middle panel: 850K, interface is rough; Bottom 
panel: 922K, interface is premelted. Al atoms are colored 
coded on the structural OP as solid (blue) and premelted liq- 
uid (green). Pb atoms are shown as red. In each image, the 
average structural OP and chemical OP are plotted as dot- 
ted green and solid white lines, respectively. The premelting 
width w is defined as the distance between the two OPs (at 
half value). 

Figure [T] shows NVT snapshots of equilibrated (111) 
Al-Pb interfaces at increasing temperature for three dif- 
ferent temperatures (625, 850 and 922K), together with 
the corresponding time-averaged structural and compo- 
sition OP profiles. For low T (t op p anel) just above T^ b , 
the (111) interface is faceted [4fJ and two profiles are 
nearly coincident. For high T (bottom panel) just be- 
low T^ 1 , the two profiles are separated by nearly 10 A 
indicating the presence of a prcmclting layer of liquid Al 
sandwiched between solid Al and liquid Pb. In the cen- 
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tral panel at 850 K, only the first complete layer of Al at 
the interface is structurally disordered. 

Results: Color-contour plots of the fine-scale density 
profiles, p(z,T), are shown in Fig. 2 for the (100) (110) 
and (111) Al-Pb solid-liquid interfaces as a function of 
temperature and distance (z) normal to the interface. 
Also, plotted in Fig.[5]are the temperature-dependent in- 
tcrfacial positions defined by the midpoints of the struc- 
tural and compositional order parameters. The distance 
between these two interfacial positions diverges as the Al 
melting point is approached, due to the formation of the 
prcmclting layer. The peaks of the density profiles are 
seen as vertical striations in the plot, which are stronger 
and more highly localized in the solid phase - to the left of 
the structural order parameter. To the right of the struc- 
tural order parameter line, the density peaks are smaller 
in magnitude and more diffuse, consistent with the usual 
structural ordering of a liquid near a surface. The much 
smaller magnitude of the liquid structural ordering in 
(110) - relative to (100) and (111) - was previously noted 
and discussed in Ref. |4(J. Note that the position of the 
interface as defined by the composition OP is roughly in- 
dependent of temperature, due to the mutual imiscibility 
of liquid Pb in both solid and premelted Al. The slight 
shift of this position towards higher z near T m for Al 
is due to the lower density (relative to the solid) of the 
growing prcmcltcd Al layer. 

Figure [2] shows that solid-liquid premelting occurs in 
this system for all orientations studied, with nearly iden- 
tical behavior. This is in contrast to surface premelting 
in many fee metal surfaces, in which (110) surfaces are 
prone to premelting, while premelting is not seen in the 
other orientations j7l.ll4|. 

Our previous simulations at 625K[40j suggested that 
the (100) and (110) interfaces are rough at that tempera- 
ture, while the (111) interface is faceted. However, visual 
inspection (together with orientational order parameter 
analysis) indicates that the (111) interface undergoes a 
roughening transition at higher temperatures (for exam- 
ple, see the middle panel of Fig. [I). To probe this behav- 
ior in further detail, we examine particle transport in the 
first layer of AL, which is characterized through the cal- 
culation of diffusion-constant profiles, D(z), determined 
from the mean-squared particle displacement versus time 
- see Ref. for details. The full diffusion constant analy- 
sis is presented in the Supplemental Information, 44 but 
here we focus on the diffusion constant measured within 
the first Al layer (that is, the particles making up the 
density peak closest to z = on the negative side in 
Fig. [2]). Figure [3] shows a log- linear plot of diffusion con- 
stant versus l/kaT for the particles in the first Al layer 
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FIG. 2. Equilibrium total number density contour maps of 
(100), (110) and (111) Al-Pb solid-liquid interface. Interface 
positions determined from structural OP profiles (circle) and 
chemical OP profiles (square) as function of T are plotted 
on top of the correspondingly contour map. The error bars 
represent 95% confidence levels. For symbols without error 
bars, the error is smaller than the size of the symbol. 



for the three interfacial orientations studied. The mag- 
nitude of the slope of this Arrhenius plot of D can be 
interpreted as an activation energy for diffusion. Except 
very close to Tj^ 1 , the slope for the rough (100) and (110) 
interfaces is constant. For (111), however, the slope un- 
dergoes a discontinuous change at a temperature of about 
826(4)K, indicating a sudden decrease in the activation 
energy for diffusion in this layer to a value that is com- 
parable to that of the rough (100) and (110) interfaces. 
This temperature where the change in activation energy 
occurs roughly where the width of the interface begins to 
show a logarithmic dependence on undercooling (see be- 
low). The temperature is also very near the roughening 
temperature of 823 K previously reported for the (111) 
interface based on in situ TEM experiments of Pb inclu- 
sions in an Al matrix. 41( The close correspondence be- 
tween these temperatures suggests that roughening and 
the onset of premelting approximately coincide for this 
interface. 



As discussed earlier, theoretical considerations 36, 
predict a logarithmic dependence of the premelting width 
on undercooling (Eq. [5]) . To examine the validity of Eq. [5] 
for the solid- liquid premelting transition in Al/Pb, we 
plot the calculated width of the premelting layer, w, as a 
function of the undercooling, AT = Tm — T on a linear- 
log plot in Fig.|4j The data is well described by Eq.[2]for 
all three interfacial orientations for undercoolings up to 




FIG. 3. Log-linear plot of the diffusion constant for Al atoms 
in the first Al layer as a function of 1/feT for the (100), (110) 
and (111) - diamonds, squares and circles, respectively. The 
Al melting point and estimated roughening temperature, T r , 
are indicated by vertical lines. The scale at the top of the 
figure shows the temperatures corresponding to the inverse 
temperatures shown on the independent axis. 



100K - deviations are seen at higher temperatures when 
the width approaches atomic dimensions, as expected. 
Using a weighted least-squares linear regression over the 
temperature range 875K to 921K, we obtain estimates for 
wo and To. Central to the derivation of Eq. [2] is the as- 
sumption that the interaction between the two interfaces 
bounding the premelted layer (the so-called "disjoining 
potential") is exponential and repulsive j4fjj|: 



4> d {w) = A 7 e- lu /™° 



(3) 



where wq is the length scale of the interaction and A7 
is given in Eq. Q] The quantity To in Eq. [2] is given by 
To = AjTm/wopL, where p is the number density and L 
is the latent heat. The fitted values of wq, Tq and A7 are 
given in Table |U This fitting gives an estimate of about 
1.4 to 1.5 A for the range of the disjoining potential that 
is relatively independent of orientation. 

Discussion and Summary: Using MD simulation we 
predict the existence of a solid-liquid interface premelt- 
ing transition at the interface between solid Al and liquid 
Pb. That is, as the melting point of Al is approached 



FIG. 4. Linear-log plot of the premelting width versus un- 
dercooling, T* 1 — T, for each of the three interfacial orien- 
tations studied. The solid lines are the results of weighted 
least-squares fits to Eq. 



TABLE I. Values of w , T and A7 (defined in Eqs. Q]and[2]) 
from a weighted linear least-squares regression of the data for 
w versus In (T m — T). Values in parentheses represent 95% 
confidence level error estimates in the last digits shown. 





(100) 


(110) 


(111) 


w (A) 
To(K) 

A7 (mJ m~ 2 ) 


1.47(3) 
1.28(12) xlO 3 
174(16) 


1.46(3) 
1.09(9) xlO 3 
148(13) 


1.36(4) 
1.45(14) xlO 3 
183(18) 



from below, the surface of the crystalline Al melts to 
form a premelting layer of liquid Al separating the solid 
Al and liquid Pb bulk phases. This transition was seen 
in the simulations for all interfacial orientations studied - 
(100), (110) and (111). Although solid- vapor and grain- 
boundary premelting transitions are well established in 
the literature, premelting of a solid-liquid interface has 
not, to our knowledge, been previously reported at am- 
bient pressures. Such a transition requires that the melt 
phase of the solid and the bulk liquid be mutually immis- 
cible, which is true for the Al/Pb system studied here. 
The width of the premelting layer is shown to depend 
logarithmically on the undercooling AT, as predicted by 
theoretical considerations. 36, III El 
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At lower temperatures, near the melting point of Pb, 
we have previously shown that the (100) and (110) inter- 
faces are rough, whereas the (111) interface is faceted - in 
agreement with experimental observations on liquid Pb 
inclusions in a solid Al matrix. (4lj In the current simula- 
tions, we observe a change in the activation energy of Al 
surface diffusion at the (111) interface at 826(4)K, which 
correlates well with both the experimental observations 
of a roughening transition at 825K for this orientation 41 
and with the observed onset of premelting in the present 
simulations (as evidenced by the logarithmic dependence 
of the premelting width on T). There are a number of 
other solid-liquid interfacial systems in which the melt 
phase of the solid is immiscible in the bulk liquid, such 
as, for example, the interface between ice and liquid hy- 
drocarbons, so further study of possible premelting in 
such systems is warranted. 
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